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We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for 
compact binary coalescence (CBC) searches during LIGO's fifth science run and Virgo's first science 
run. We present noise spectral density curves for each of the four detectors that operated during 
these science runs which are representative of the typical performance achieved by the detectors 
for CBC searches. These spectra are intended for release to the public as a summary of detector 
performance for CBC searches during these science runs. 



I. INTRODUCTION 



In this note, we summarize the sensitivity to gravitational waves from compact binary coalescences (CBCs) achieved 
during LIGO's fifth science run (S5) and Virgo's first science run (VSRl) [HE]- A complete and search-independent 
characterization of the sensitivity of a gravitational wave detector over short timescales is given by the spectral 
density. Over long timescales, the noise in LIGO and Virgo detectors is non-stationary and a spectral density is not 
an appropriate description of a detector's sensitivity. Our goal here is to characterize the overall performance of each 
detector for CBC searches in all of S5 and VSRl using the familiar notion of spectral density. However, since the 
concept of spectral density is not appropriate for LIGO and Virgo detectors on such a long timescale, we choose a 
single "representative" time for each detector and compute the spectral density in each detector around this time. 
We take the resulting spectral density as representative of the typical sensitivity to gravitational waves from CBCs 
achieved by the detectors in S5/VSR1. 

We choose the representative spectral density curves to correspond to times when each detector operated near the 
mode of its inspiral horizon distance distribution. The inspiral horizon distance is a quantity derived from the spectral 
density that summarizes the sensitivity of a detector to gravitational waves from CBCs at a given time. In this article, 
we gather the inspiral horizon distance data generated during S5/VSR1 CBC analysis and use the results to identify 
spectral density curves that are representative of detector performance for CBC searches during these science runs. 

The plots and data presented here are intended to be released to the public as a summary of detector performance 
for CBC searches during S5 and VSRl. These results use exactly the same science segments and analysis code [3] 
that was used in the CBC searches in S5 and VSRl. The results presented here supersede previous studies g] done 
on the inspiral horizon distance in that we use LIGO's version 3 calibration [S], the same calibration version used in 
the S5 CBC search. Using version 4 calibration |S] would change these results by about 5 to 10%, but not in a simple 
way. Since version 4 h{t) has never been made for S5 and version 3 calibration was used in the S5 search, we stick 
with version 3 calibration here. However, we do include first order corrections for version 4 DC calibration, which 
affects the results through an overall multiplicative scaling factor. All CBC searches applied these DC calibration 
corrections, so the results presented here compare directly with the CBC search sensitivities reported in our papers. 
For Virgo, we use version 2 calibration j7j , the same calibration used in the VSRl CBC search. 
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In the next section, we define the inspiral horizon distance and present a summary of the inspiral horizon distance 
data from S5/VSR1 CBC analyses. In section 3, we explain how we used the inspiral horizon distance data to compute 
a spectral density that is representative of each detector's sensitivity for CBC searches in S5/VSR1. 



II. INSPIRAL HORIZON DISTANCE 



The (power) spectral density S'„(/) for a stationary random process n(t) is defined implicitly by the relation 

lsMm~f)^{n{f)n*{n), (1) 

where n{f) is the Fourier transform of the random process. The spectral density is a measure of the mean square 
noise fluctuations at a given frequency. In LIGO and Virgo applications, we treat the strain noise in a detector as a 
stationary random process. If the noise in the detector were truly stationary, then the noise spectral density would 
completely characterize the sensitivity of the detector as a function of frequency. 

As mentioned above, the noise in the LIGO and Virgo detectors is not stationary. However, by measuring the 
spectral density over a short enough timcscale, we are able to approximate the noise as stationary. The chosen 
timescale must also be long enough that we can form an accurate estimate of the spectral density. In the S5/VSR1 
CBC searches, the spectral density was computed on 2048-second blocks of contiguous data [5j. We account for long 
timescale non-stationarities by using a different spectral density for every 2048 seconds. 

In assessing the overall performance of a detector for CBC searches, we use the inspiral horizon distance data from 
S5 and VSRl to identify the "typical" sensitivity of the interferometers. The inspiral horizon distance of a detector 
is the distance at which an optimally oriented and optimally located equal-mass compact binary inspiral would give 
an average signal to noise ratio (SNR) of p = 8 in the interferometer. If h{f) represents the Fourier transform of the 
expected signal, then the average SNR this signal would attain in a detector with spectral density Sn{f) is given by 




(.) = >/4/ If./. (2) 

We find the inspiral horizon distance by setting (p) = 8 and solving for the distance D to the inspiral event which 
parametrizes the waveform h(f). Thus, the inspiral horizon distance combines the spectral density curve with the 
expected inspiral waveform to produce a single quantity that summarizes the sensitivity of the detector at a given 
time. 

Practical considerations require modifications to the limits of the integral. In the CBC search code, we compute 
the signal to noise ratio by 




The lower limit is determined by our ability to characterize the noise at low frequencies. In the S5 CBC search, we 
took flow = 40IIz as the low frequency cut-off in computing the inspiral horizon distance. For Virgo in VSRl, the 
low frequency cut-off was flow — 60IIz. The upper limit of the integral is the innermost stable circular orbit (ISCO) 
frequency, 

^'^'"^ 6V67TGM' ^^'^ 

where M is the total mass of the binary system. For binary neutron star systems, fisco — 1570Hz. However, the 
inspiral horizon distances reported here are culled from the S5 high mass search (except for Virgo), which down- 
sampled the h{t) data to 2048 Hz, so that the integral was cut off at f^y — 1024 Hz. However, this is a < 1% effect, 
even in computing the inspiral horizon distance for low mass systems, because most of the SNR is accumulated in the 
"bucket" of the noise curve; we neglect it here. 

The inspiral waveform for CBCs is accurately given in the frequency domain by the stationary phase approximation. 
For an optimally oriented and optimally located equal mass binary, the signal that appears at the interferometer (in 
this approximation) is given by 

Hf) = ^ (^) (GA^)^/«(7r/)-V6e»*(/;M)^ (5) 
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Inspiral Horizon Distance During S5 and VSRl 
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FIG. 1: Distribution of inspiral horizon distance for the four gravitational wave detectors HI, LI, H2 and VI for all of S5 
and VSRl. This histogram includes each 2048-second analyzed segment from S5 and VSRl. The distributions shown here 
correspond to the 1.4 -1.4 solar mass inspiral horizon distance for the LIGO detectors. For the Virgo detector, we have plotted 
the 1.0-1.0 solar mass inspiral horizon distance distribution, scaled by (2.8/2)^''^ to adjust for the lower mass. 



where M is the chirp mass of the binary, D is the distance to the binary and is a real function of /, parametrized 
by the total mass M. Setting (p) = 8 and inserting this waveform into eqn. [sj we find that the inspiral horizon 
distance is given by 



8 I 24c3 / 



1/2 



(GM) 



5/6^-7/6 




gh f-7/3 



df, 



(6) 



where D is expressed in Mpc. The inspiral horizon distance is defined for optimally located and oriented sources. 
For a uniform distribution of source sky locations and orientations, we divide the inspiral horizon distance by 2.26 to 
obtain the SenseMon range [9] reported as a figure of merit in the LIGO and Virgo control rooms. 

In practice, it is convenient to measure distances in Mpc and mass in Mq. It is useful therefore to specialize eqn. 
[6] to this unit system. Further, since we measure the strain h(t) at discrete time intervals At = l//s, the spectral 
density is only known with a frequency resolution of A/ = fs/N, where N is the number of data points used to 
measure 5'„(/). By putting / — kAt into eqn. |6]and grouping terms by units, we arrive at the expression 



\ 



4 
N 



(7) 



where 



r 



\ 247r4/3 



) ^v^j [m^J [ IMpc J [ At 



-1/6 



(8) 



for the inspiral horizon distance in Mpc. Since it is convenient to work with the binary system's component masses, 
we have also replaced the chirp mass M with the reduced mass fi and the total mass M, where M = p^/^M"^/^. 
Written this way, the inspiral horizon distance in Mpc is easily computed from the binary component masses in Mq. 
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Inspiral Horizon Distance During S5 and VSRl 
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FIG. 2: Mean inspiral horizon distance as a function of mass for the four gravitational wave detectors HI, LI, H2 and VI 
during all of S5 and VSRl. The error bars on the curves extend from one standard deviation below to one standard deviation 
above the mean. 
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FIG. 3: Inspiral horizon distance as a function of time during S5 and VSRl. The average inspiral horizon distances for each 
week in S5 and VSRl. As an indication of the weekly variations, we have included error bars corresponding to the standard 
deviation of the inspiral horizon distance during each week. 
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TABLE I: Parameters used in the computation of the spectral density. 



FINDCHIRP parameter 'W 


S5 high-mass 


VSRl low-mass 


representative spectra 


sample rate (l/At) 


2048 Hz 


4096 Hz 


16384 Hz 


data block duration (Tuock) 


2048s 


2048s 


2048s 


number of data segments (Ns) 


15 


15 


1023 


data segment duration (T) 


256s 


256s 


4s 


stride (A) 


262144 


524288 


32768 



We have collected the inspiral horizon distance data from the original output of the S5 high-mass and S5/VSR1 
low-mass CBC searches to be published. We have collected the data, rather than computing the inspiral horizon 
distance directly, in order to ensure that we analyze the exact same science segments and use the exact same analysis 
code as used in the LIGO/Virgo CBC searches. 

In fig. [l] we histogram the inspiral horizon distance for the four gravitational wave detectors HI, LI, H2 and VI 
for all of S5 and VSRl. The plot depicts the variability in sensitivity to gravitational waves from CBCs throughout 
the science runs. The bimodal behavior seen in the LIGO detector distributions is due to a significant commissioning 
break in S5, which greatly improved the sensitivity of the LIGO detectors. 

In fig. [2] we plot the mean inspiral horizon distance for each interferometer as a function of the binary component 
mass. This plot reflects the mean performance of the detector over various frequency bands. As the component 
mass becomes higher, the upper cutoff frequency fhigh — fisco becomes smaller and smaller. This means that the 
inspiral horizon distance focuses on a narrower band around the lower cutoff //o^ — 40Hz (or flow — 60IIz in the case 
of Virgo). The fall-off of the inspiral horizon distance for high-mass binaries characterizes the performance of the 
detector near the low cut-off low frequency. In particular, it should be emphasized that the sensitivity of the detectors 
during S5 and VSRl to high-mass systems does not fall off as might be suggested by the graph. The inspiral horizon 
distance takes into account only the inspiral stage of the CBC event, while for high-mass systems (M > 25Msun) the 
merger and ringdown stages of the occur in the LIGO and Virgo sensitive band. For these binary systems, we use 
EOBNR waveform templates that include the merger and ringdown stages and our sensitivity is significantly greater 
than depicted here. 

For purely historical reasons, the S5 high-mass runs, which did not include VI, computed the inspiral horizon 
distance for (n-f 0.4)-(7i-|-0.4) solar mass binaries for integers n > 1. The S5 low-mass runs, which included VI, 
instead computed the horizon for n-n solar mass binaries. In order to make an apple-to-apples comparison, we scaled 
the Virgo distribution by (2.8/2)^/^ corresponding to the ratio of chirp masses for the LIGO and Virgo data. This 
scaling ignores the fact that fisco is different for the two mass pairs, but this is negligible since the template is buried 
in the noise at that high of a frequency. 

III. REPRESENTATIVE NOISE SPECTRAL DENSITY 

In this section, we present spectral densities which we suggest are representative of the sensitivity achieved for 
S5/VSR1 CBC searches. The chosen representative curve corresponds to a time when the detector operated near the 
mode of its inspiral horizon distance distribution. We also illustrate the variability in detector sensitivity by giving 
spectra for HI corresponding to times when HI operated near the mode, early-S5 mode, mean and max of its inspiral 
horizon distance distribution. 

The algorithm used to compute the spectral densities is described in detail in [S]. The parameters needed in order 
to reconstruct our results and are given in table |l] The first column in table [l] gives a list of parameter names and 
symbols, which are the same names and symbols used in [8]. The second and third columns gives the values of 
these parameters used in S5/VSR1 CBC searches. These parameters can be used to reproduce the inspiral horizon 
distance data accompanying this note. The fourth column gives the values of the parameters used to compute the 
representative spectral density curves shown here. In making our choice of parameters for computing representative 
spectra, we sacrificed frequency resolution (A/ — l/T) for PSD accuracy (which increases with Ng). 

In fig. |4] we give representative spectral density curves for each of the four detectors during S5 and VSRl. These 
curves correspond to times when the detectors operated near the mode of their 1.4-1.4 solar mass inspiral horizon 
distance distribution. The strain calibration is valid for LIGO above 30 Hz and for Virgo above 10 Hz, which give the 
lower limits for the curves plotted here. 

In table|Tlj we provide a quantitative summary of the low-mass inspiral horizon distance distributions. All distances 
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Representative Noise Spectral Density Curves During S5 and VSRl 

10-13 , . I ■ ' I ■ ■ — ^ 




10"" I ■ ■ — • ' — ' 

10^ 10= 10^ 

Frequency (Hz) 



FIG. 4: Representative spectral density curves for LIGO and Virgo detectors during S5 and VSRl. These spectral density 
curves correspond to July 21, 2007 (GPS 869040003) for HI, March 16, 2007 (GPS 858087370) for LI, September 30, 2007 
(GPS 875211248) for H2 and June 13, 2007 (GPS 865748914) for VI. These times are chosen such that the inspiral horizon 
distance for each detector at that time coincides with the mode of the its inspiral horizon distance distribution, as given by the 
midpoint of the most populated bin in fig. [l] 



TABLE IL Summary of Inspiral Horizon Distance Data from S5 and VSRl 





HI 


LI 


H2 


VI 


mean 


30.4 


30.3 


14.1 


7.8 


majc 


36.2 


37.5 


16.7 9.8 


mode 


34.1 


33.1 


15.3 


7.7 


std 


3.5 


4.5 


1.9 


0.6 



are given in Mpc. The Virgo data was for 1.0-1.0 solar mass but has again been rescaled by (2.8/2)^/^. These statistics 
give varied viewpoints on the detector performance. In the case of LI, for instance, the mode and max differ by more 
than 10%. This fact suggests that the spectral density curves corresponding times when LI operators at its mode 
differ significantly from its "best" spectral density. 

To illustrate this point, we plot in fig. [5] four spectra for HI from different times in S5. We see a significant 
improvement in the spectral density from the beginning of S5 to the end of S5. In early S5, HI operated most often 
with a 1.4-1.4 solar mass inspiral horizon distance near 26.7 Mpc. In late S5, the inspiral horizon distance distribution 
peaked around 36.6 Mpc (early S5 roughly corresponds to the lower peaks in the horizon distributions from fig. [T] 
while late S5 roughly corresponds to the higher peaks). Correspondingly, we see a decrease in the spectral density 
from early S5 to late S5. 

All of the LIGO data used here has been computed using version 3 calibration. Updating these results using version 
4 calibration would be nontrivial. However, we have corrected all plots and data presented here for changes in the DC 
calibration from version 3 to version 4, as reported [Hj. This correction amounts to multiplying the spectra for each 
detector by an overall scaling factor. The scaling factors for the various detectors are 1.07 for HI, 1.062 for H2 and 
0.96 for LI. The inspiral ranges are reduced by this same factor. These corrections make HI and H2 are somewhat 
less sensitive in version 4 calibration, while LI is somewhat more sensitive. Note that the noise spectra presented 
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Representative Noise Spectral Density Curves For HI During S5 
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FIG. 5: Representative spectral density curves for HI during S5. These spectral density curves correspond to times when 
the detector operated near its late S5 mode (34.1Mpc), early S5 mode (24.8Mpc), mean (30.5Mpc) and maximum (36.2Mpc) 
inspiral horizon distance. The times chosen are July 21, 2007 (GPS 869040003) for the late S5 mode, February 05, 2005 (GPS 
823205705) for the early S5 mode, January 25, 2007 (GPS 853767368) for the mean and August 15, 2007 (GPS 871198828) for 
the maximum. 



here are subject to systematic uncertainties associated with the strain calibration. These uncertainties can be up to 
±15% in amphtude. For more detail, see references [Sl[7]- 



IV. CONCLUSIONS 



We have presented spectral density curves for each of the four detectors used in S5 and VSRl. We suggest that 
these noise curves are representative of the sensitivity of the LIGO and Virgo detectors for CBC searches in S5 and 
VSRl. Our choice for these noise curves is based on the inspiral horizon distance, which is a measure of detector 
performance for binary inspirals. We intend for these noise curves to be released to the public as a summary of 
detector performance during S5 and VSRl. 
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